First-principles-based ±s-wave modelling for iron-based superconductors: Studies for 
specific heat and nuclear magnetic relaxation rate 
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In order to consistently explain controversial experimental results on superconducting states ob- 
served by diflterent probes in typical iron-based superconductors, we construct a realistic multi-band 
±s-wave pairing model by combining the quasiclassical formalism with the first-principles calcula- 
tion. The model successfully resolves the controversies in contrast to the fact that simplified models 
such as two-band ±s-wave one fail to do. A key in the model is the existence of relatively small 
gaps which leads to material-dependent peculiarities. 

PACS numbers: 74.20.-z, 74.70.Dd, 74.25. Jb 



The discovery of the iron-pnictide superconductor 
LaFeAs(Oi_:,F:,) [H made a striking impact on materi- 
als science, because this compound includes the element 
of the most familiar ferromagnetic metal, Fe, as a main 
component. The transition temperature in the so- 
called "1111" compounds i?FeAs(Oi_:rFx) (i?==Pr, Nd, 
Sm) exceeds 50K, which is the highest except for high 
cuprates. In addition to the high transition temperature, 
the variety of related materials is quite rich. For exam- 
ple, "122" compounds (Ai_a;B^)Fe2As2 (A=Ba, Sr, Ca, 
B=K, Cs, Na) and "11" compounds Fe(Se2;Tei_x) are 
the typical family materials, whose element substitutions 
are widely possible @, particular, the supercon- 

ductivity is surprisingly robust against substitutions of 
Ni and Co for Fe. 

In contrast to the discovery rush of family compounds, 
their superconducting states still remain elusive. There 
is no established pairing-symmetry model explaining all 
experimental results consistently. In the early days, puz- 
zling experimental data were reported. In spite of sharp 
resistivity drops and clear Meissner signals at the super- 
conducting transition Tc, the jump of C/T, where C is 
the specific heat and T is the temperature, was hardly 
observable in 1111 compounds. The reason was initially 
ascribed to a large phonon contribution which masks an 
electronic one. Afterwards, a small jump at Tc and the 
concave-down temperature dependence of C /T below Tc 
were confirmed in not only 1111 compounds 0, Q but 
also structurally equivalent LaFeP(0i_2;F2:) [6[. This 
fact strongly suggests the existence of rather small su- 
perconducting gap Q in addition to large main gaps in 
the cases of s-wave gap. In this paper, we clarify that 
such a multi-gap structure consistently explains all ex- 
perimental observations of 1111 compounds by means of 
a realistic five-band is-wave pairing model based on a 
first-principles calculation. A striking result of the model 
is a natural reproduction of the nuclear magnetic relax- 



ation rate 1/Ti below Tc- 

On the other hand, 122 compounds experimentally ex- 
hibit large jumps and exponential behavior in C/T like 
conventional superconductors @, [l,S[l3|- Moreover, the 
power law in the T-dependence of below Tc is differ- 
ent from that of 1111 compounds [1, [ill, 12,[l3|- In fact, 
an angle resolved photoemission spectroscopy (ARPES) 
study reported that all gaps fully open and the differ- 
ence between their gap amplitudes is not so significant 
[l^ . In this case, we reveal that the multi-gap structure 
according to ARPES data can consistently explain the 
specific heat and 1/Ti data without any assumptions ex- 
cept for ±s-wave paring symmetry [isl . [l6l |. From the 
present analyses on 1111 and 122 compounds, it is found 
that the existence of the relatively small gap gives rise to 
the material variety. We believe that this fact has an key 
role on the quest for the superconducting mechanism. 

Let us present a procedure to construct the realistic 
multi-band model. For the band structure around the 
Fermi level Ep, we perform a first-principles calculation 
[l7| . which provides multiple Fermi-surfaces and their 
density of states (DOS) at Ef depending on the tar- 
get materials. When evaluating superconducting gaps 
for multi-band superconductors, we examine all the ex- 
perimental data and select key data. ARPES measure- 
ments are successful for 122 compounds [l3|- Therefore, 
we can directly fit ARPES data to determine the gap 
amplitudes on each band in the 122 case. On the other 
hand, ARPES [3| measurements are technically difficult 
for 1111 compounds. Instead, to estimate the gap am- 
plitudes we adopt data oi C/T Q and the penetration 
depth 3, 2^ in the 1111 case. In fact, those data clearly 
suggest that the pairing symmetry is full-gap but a small 
sing le or small multi gaps coexist with main large gaps 
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2^,l23|,l24|. 



We describe the expression for C/T here, while we refer 
readers to Refs. (HSH^I for 1/Ti. The T-dependence 
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FIG. 1: (color online) (a) The band structure calculated by the generalized gradient approximation using structural measure- 
ment values of BaFe2As2. (b) The Fermi surfaces and the density of states at the Fermi energy. Indices {a—1, 2, 3, 4, 5) are 
assigned from r(zone center) to X. Temperature dependences of the superconducting pair-potential are displayed in (c) 
and (d). Here, "B" of Case Bl or B2 stands for BaFe2As2. (e) Temperature dependences of (Cs — Cn)/T. Cs(n) is the specific 
heat of the superconducting state (normal state), (f) Temperature dependences of the nuclear magnetic relaxation rate 1/Ti. 



of C/T is calculated by the second derivative of the free 
energy. The free energy can be evaluated by the qua- 
siclassical theory of superconductivity [H, [2g| , which is 
a mean field treatment convenient in evaluating super- 
conducting properties. To calculate C/T and 1/Ti, we 
need the T-dependence of the multiple superconducting 
pair-potential on each band. For this purpose, we 
solve the gap equations for multi-band superconductors 
[sot . With Matsubara frequency cOn = (2n -|- 1)ttT, the 
gap equations are written as 



Aa = 27rT ^ Kpfpii 

w„>0 /3 



(1) 

(2) 



where a and /3 stand for the band index, and Xap is the ef- 
fective coupling constant. In addition, Xap — NpXpa/Na 
with Na being DOS at for a-band. Xaa comes from 
the intra-band interaction, and Xap, where a ^ P, gives 
the pair- hopping between the different bands. The effec- 
tive coupling constants Aq/j's are employed as parameters 
to reproduce the experimental results for the supercon- 
ducting gaps. Again, we note the selected experiment 
typ e depending on a kind of the compounds, i.e., ARPES 
[1J| for the 122 compounds, and the specific heat ^6] and 



the penetration depth 



for the 1111 compounds. 



and f^{iujn) follow the Eilenberger equations as, 

UJnfa{il^n) = Aaga{il^n), (3) 

Wn/I(iw„) = A^gaiiuin), (4) 

giiiuJn) = 1 - fa{iuJn)fl{i(^n), (5) 

where Re ga{iujn) > for a;„ > 0. The free energy dif- 
ference between the superconducting and normal states 

iillillil,!!!, F,n = F,, 



super -^normal; GXpiGSSCd &S 

1 - 5Q(iw„) 



1 -I- ga{itOn) 



A*/a(lW„) 



which demands the solutions of Eqs. (3), (4), and (5) 

ga{iLOn) 



(6) 

(7) 

(8) 
(9) 



Note that /^^•'(iLJn) and gQ,(ia;„) — > 1 when 
ujn ^ |Aq|. The cut-off frequency u^is introduced as 
X;";>o in Eqs. (1), (2) and (6) 0. Equation (6) 
indicates that Fgn can be directly evaluated by DOS Na 
obtained from first-principle calculations and the gap 
values Aq,. 



The quasiclassical Green's functions ga{i(-Un), faiii-^n 
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FIG. 2: (color online) (a) The band structure calculated by the generalized gradient approximation using structural measure- 
ment values of LaFeAsO. (b) The Fermi surfaces and the density of states at the Fermi energy. Indices (q=1, 2, 3, 4, 5) are 
assigned from r(zone center) to M. Temperature dependence of the superconducting pair-potential Ac in (c) and (d). Here, 
"L" of Case LI or L2 stands for LaFeAsO. (e) Temperature dependences of (Cs — C-n)/T. Cs(n) is the specific heat of the 
superconducting state (normal state), (f) Temperature dependences of the nuclear magnetic relaxation rate 1/Ti. 



Since C/T at constant volume is generally obtained 
by C/T = dS/dT = -d^F/dT^, the specific heat in the 
superconducting state is expressed as 



Cs — C„ 
T 



QJ'2 



(10) 



where Cs{n) is the specific heat of the superconducting 
state (normal state). We can rewrite C^/T as 7n, which is 
the Sommerfeld coefficient. Thus, we numericall y c alcu- 
late d^Fsn/dT^ to obtain Cs/T. As for l/Ti [H, Hi, [13] , 
we calculate it using Na and Aq, obtained here. 

Let us present the calculated results. The first focus 
is 122 compounds. Figure 1(a) shows the band struc- 
ture of BaFe2As2. This result is obtained by using the 
generalized gradient approximation (GGA) based on the 
measured structural data Q. The obtained Fermi sur- 
faces are displayed in Fig. 1(b). The calculation gives 
DOS's at Ep, which are input parameters in the gap 
equations (1) and (2) and the free energy (6). Each 
DOS at Ep is written as Na, where a is numbered as 
a = 1, 2, 3, 4, 5 from F point (zone center) to X as shown 
in Fig. 1(b), and a = 1 to 3 (4 to 5) correspond to hole 
(electron) bands. Throughout this paper, we adopt ±s- 
wave pairing model since other choices fail to reproduce 
the experimental data consistently. The positive (nega- 
tive) sign is assigned to Aq, of hole (electron) bands (Aq, 
are assumed to be real). The T-dependence of Aq is cal- 
culated so as to follow the ARPES result Q (Case Bl) 
as shown in Fig. 1(c). On the other hand, in Fig. 1(d) 



we show, for comparison, the result (Case B2) with the 
minimum jAaj being slightly bigger than the ARPES 
result. Figure 1(e) shows the T-dependences of C/T. 
Both Cases Bl and B2 show no significant difference from 
the weak coupling single-band Bardeen-Cooper-Schrieffer 
(BCS) result. One of the reasons is that the weighting 
of the small-gap band is small compared to total one, 
i.e., the ratio Ns/Nt is 0.2, where Nt is the total DOS 
{Nt = J2a^a)- Moreover, the gap-amplitude difference 
between the minimum jAsj and the maximum IA2I is 
not so large, where IA3I/IA2I - 0.5 (-- 0.7) for Case Bl 
(B2). On the other hand, the influence of |A2|/Tc ^ 2, 
which is bigger than the single-band's value 1.76, slightly 
enhances C/T at Tc. As a result, the difference in C/T 
becomes small among Case Bl, B2 and single-band case. 
The T-dependences of C/T in Cases Bl and B2 are con- 
sistent with the experimental observations 0]. 

On T-dependence of 1/Ti, both the cases show signifi- 
cant differences from the single-band BCS case as shown 
in Fig. 1(f). The coherence peak just below Tc is absent 
in both Cases Bl and B2, because the cancellation be- 
tween and "— " signs of Aq, is effective. In contrast, 
in the single-band BCS case, even if the damping rate of 
the quasiparticle is taken large as 77 = O.lTc, the peak is 
clearly identified as seen in Fig. 1(f). Moreover, the low- 
lying excitation arising from the small gap IA3I alters T- 
dependence of 1/Ti compared to the single-band s-wave 
case. We point out that the five-band model Case Bl 
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successfully reproduces the experimental results of C/T 
and 1/Ti 

We next turn to 1111 compounds. At first, based on 
GGA with measured structural parameters [3J|, we ob- 
tain the band structure of LaFeAsO as shown in Fig. 2(a). 
The band index is also numbered as a = 1, 2, 3, 4, 5 from 
r (zone center) to M. The five Fermi surfaces are dis- 
played in Fig. 2(b), where TVi, N2, N3 (N^, N5) are the 
DOS at Ep of hole (electron) bands. We prepare two 
types of multi-gap structures, Cases LI and L2, whose 
T-dependences of are plotted in Figs. 2(c) and 2(d), 
respectively. These are estimated by experimental data 
of C/T 0,1, Q and the penetration depth [i^], because 
direct experimental ARPES data of the gaps are not 
presently available in 1111 compounds. Case LI gives 
four large gaps and one small gap, while Case L2 consid- 
ers a medium gap between the maximum and minimum 
gap-amplitudes. Figure 2(e) shows T-dependence of C/T 
for both Cases LI and L2. When the amplitude differ- 
ence between the minimum and maximum gaps is large, 
the jump of C/T at T^ decreases. For comparison, we 
note that the small gap in the 122 compounds as shown 
in Fig. 1(e) is not enough to reduce the jump of C/T 
like the present 1111 compound case. In addition. Case 
L2 reproduces the concave-down behavior as observed in 
LaFeP(Oi_a,Fa;) Thus, the small-jump feature and 
the concave-down behavior in C/T suggest not only the 
existence of a small-gap but also that of a medium-gap 
band, whose contributions are significant compared to 
122 compounds. We then expect that Case L2 is the 
most possible candidate. 

Figure 2(f) shows T-dependence of 1/Ti for both Cases 
LI and L2. We find that both Cases LI and L2 surpris- 
ingly exhibit T^-behavior of 1/Ti up to the experimen- 
tally accessible low T. The low-lying excitations due to 
the small gap push the exponential behavior of 1 /Ti into 
a further lower temperature region. Moreover, the peak 
just below Tc does not appear because of the cancellation 
due to ± signs The most of 1111 compounds show 
T3 dependence in 1/Ti below Tc [1113], which suggests 
that Case L2 is the best as noted in Fig. 2(f). 

In conclusion, we examined the validity of ±s-wave 
scenario for typical iron-based superconductors (122 and 
1111 compounds) through the realistic model using the 
quasiclassical formalism combined with first-principles 
calculations. Consequently, we found that any anoma- 
lous properties observed in the specific heat and the 
nuclear magnetic relaxation rate are reproducible with- 
out any extrinsic assumptions, i.e., all required is the 
properly-evaluated gap amplitude of each band. 

Finally, we add a note that the momentum dependence 
of the gap amplitudes is reported by some recent ex- 
periments [H m, [111 of the related compound. Such 



an anisotropy can be easily implemented in the present 
framework. 
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